Characterizing optical chirality 
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We examine the recently introduced measure of chirality of a monochromatic optical field [Y. Tang 
and A. E. Cohen, Phys. Rev. Lett. 104, 163901 (2010)] using the momentum (plane-wave) repre- 
sentation and helicity basis. Our analysis clarifies the physical meaning of the measure of chirality 
and unveils its close relation to the polarization helicity, spin angular momentum, energy density, 
and Poynting energy flow. We derive the operators of the optical chirality and of the corresponding 
chiral momentum, which acquire remarkably simple forms in the helicity representation. 
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I. INTRODUCTION 

Optical chirality plays an important role in the optical 
sensing of biomolecules (see, e.g., [1]) and the interac- 
tion of light with chiral nanostructures or metamateri- 
als [2]. The nonreciprocity in transmission or circular 
dichroism measures the chiral effects in the interaction 
of the optical field with the specimen. However, it is of 
both fundamental and practical interest to introduce a 
measure of the chirality of the optical field itself. Re- 
cently, Tang and Cohen [3] used a proposal by Lipkin 
[4] and introduced the local measure of the chirality of 
a nonparaxial monochromatic field. This will be called 
here optical chirality density in order to emphasize its 
local nature. Together with the chirality flow density, 
these quantities satisfy the continuity equation, akin to 
the Poynting theorem. A recent experiment [5] confirmed 
that this measure of optical chirality is meaningful. 

One of the points of [4] was that the chirality density is 
not simply related to the local ellipticity (i.e., the degree 
of circular polarization of light in real space) but repre- 
sents a more sophisticated characteristic which can take 
arbitrary large values when divided by the electric field 
energy density. Here we show that the chirality density 
is directly related to the helicity of light, i.e., the degree 
of circular polarization in the momentum (plane-wave) 
representation. We derive local and integral values of 
the chirality and chirality flow of a nonparaxial free field 
and show that the operators of energy and momentum 
of light, multiplied by the helicity, correspond to these 
quantities. This unveils the actual relation of the opti- 
cal chirality to the polarization and sheds light on the 
similarity with the Poynting theorem. 



sian units, we write the energy density w and Poynting 
energy flow s [6]: 



w = 9 -(£ 2 + n 2 ) 



s = cg(£ x U) 



(1) 



(2) 



where g = (47r) 1 . The chirality density \ an( i the cor- 
responding chirality flow ip introduced in [3,4] read: 

x = - § [s ■ v x e + n • v x u\ , (3) 
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if = — Z [£ x (V x U) - U x (V x £)) . (4) 

uj 2 

Compared to [3], we multiplied Eqs. (3) and (4) by a 
constant c/uj, in order to have the same dimensionality 
as the energy density and flow. The energy and chirality 
satisfy the continuity equations [3,6]: 



dw 
~dt 

dx 
dt ' 



V-s = , 



V • (p = 



(5) 
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It is worth noticing that \ an d V 3 are time-independent 
[4] (so that dx/dt = 0), whereas w and s possess oscil- 
lating terms [6]. Performing time averaging, we express 
these quantities via complex fields: 



w = - t Re (E* • E + H* • H) 
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(7) 



II. BASIC EQUATIONS 

To begin with, we consider monochromatic light in 
free space, characterized by the real electric and mag- 
netic fields: £(r,t), H(r,i), and their standard com- 
plex representations: E(r), H(r), so that £(r,t) = 
Re [E (r) e~ iut ] , U{r, t) = Re [H (r) e" 1 "*] . Using Gaus- 



s = c I Re (E* x H) , (8) 



£ = X = -|lm(E*-H) , (9) 



ip = if = c | Im (E* x E + H* x H) . (10) 
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Here we used Maxwell equations and wrote Eqs. (7) 
(10) in a form exhibiting a notable symmetry between 
the energy and chirality. While w and s are respectively 
a scalar and a vector, \ an d ¥> are a pseudoscalar and 
a pseudovectors changing their signs upon mirror reflec- 
tions. After performing the time-averaging, the continu- 
ity equations (5) and (6) reduce to V-s = and V-<p = 0. 



III. HELICITY REPRESENTATION 

It is known that the chirality of an electromagnetic 
field is typically associated with the degree of circular 
polarization. At the same time, this quantity, known 
as helicity (i.e., the spin state of a photon), is well de- 
fined only for plane waves, i.e., in the momentum rep- 
resentation [7]. This is due to the transverse nature of 
the electromagnetic waves, for which the polarization is 
orthogonal to its wave vector. Therefore, one can nat- 
urally characterize the polarization of each plane wave 
in the spectrum, but it is difficult to characterize the 
local spatial polarization of a nonparaxial superposition 
of multiple plane waves - in the generic case, the po- 
larization has all three components and exhibits rather 
complicated features [8]. Recently, analyzing the energy 
flows in optical fields, we have shown that the dynamical 
characteristics (such as energy, momentum, and angu- 
lar momentum) of a generic propagating electromagnetic 
field acquire a particularly clear and simple form in the 
helicity momentum representation [9]. Below we apply 
this representation to Eqs. (7)-(10). 

The Fourier plane- wave spectrum of the complex fields 
can be written as 



{E(r),H(r)} = |- J |e (k) , H (k) j e ikr d 2 k , (11) 

where the integration is performed over the (kxiky) 
plane (for simplicity we assume propagating field and 
neglect evanescent waves) and the normalization factor 
a = \j2oj~fg is introduced for convenience below. Now 
each complex Fourier amplitude can be represented as 
a sum of two circularly-polarized plane waves with well- 
defined hclicities: 



E = E+ + E , H = H+ + H 



(12) 



where the helicities a = ±1 correspond to the right-hand 
and left-hand circularly-polarized waves. The basic prop- 
erties of the circular polarizations yield [7] 



H a = -iaYf , E CT * x E" = if 



(13) 



Substituting the representation (11)— (13) and perform- 
ing some calculations similar to those in [9,10] for the 
energy flows, we arrive at 



w 



=±i 



d 2 k' I d l ke 1 ' 



(e<"' • E ff ) 



^2 Im f d 2 k' f d 2 ke iK T a (e ct *' x E ff ) , 
(2tt) a=±1 J J 

(15) 



X = V Re / d 2 k' f d 2 ke iKr a (e ct *' • E ff ) 

( 27r ) a =±l J J 



(16) 



cp = 



CUJ 



Im J d 2 k' J d 2 ke lKr (E CT *' x E CT ) , 



(17) 

where k = k - k' and E CT *' = E'™ (k'). It is clear that in 
Eqs. (14)-(17) the integrands of the energy and chirality 
quantities differ only by the helicity a, which flips upon 
mirror reflections and spatial inversion. Note also that 
Eqs. (14)-(17) do not contain interference cross-terms 
mixing different helicities - this is a remarkable feature 
of the helicity representation diagonalizing quadratic field 
forms of Maxwell equations [9,10]. 

From Eqs. (14) and (16) it immediately follows that 
the ratio of the chirality density \ to the energy density 
id cannot exceed 1 in absolute value: 



X t» T — w ' 

W W + + W ' 



-1,1] • 



(18) 



(14) 



Here w' 7 is the energy density of the er-polarized part of 
the field and the maximal chirality x/w = c is achieved 
for a field composed of plane waves with the same helicity 
a. It should be emphasized that in this case the actual 
local polarization of the optical field in real space (which 
results from interference of multiple plane waves propa- 
gating in different directions) can be rather complicated 
and far from being circular [8,10]. 

Our results clarify the examples given in the Support- 
ing Online Materials to [3], where a superposition of 
two counter-propagating plane waves with different po- 
larization was considered. One can note that the chiral- 
ity density obtained for such field represents the sum of 
energy-weighted helicities of the partial waves. The enor- 
mously high response in the so-called superchiral fields 
was achieved in [3,5] , because there the chirality efficiency 
was determined by the ratio of the chirality density to the 
electric field energy density, which is only a part of the to- 
tal energy density of the field. (Indeed, in the example [3] 
of counterpropagating orthogonally-polarized waves, the 
maxima of the electric energy density correspond to the 
minima of the magnetic energy density and vise versa.) 
Note also that in [5] the evanescent near-field modes were 
involved, which are beyond the scope of this work. 

It is interesting to calculate the integral values of the 
energy, momentum, chirality, and corresponding 'chiral 
momentum' of the field. The momentum density p differs 
from the energy flow s by a factor of c 2 (p = s/c 2 ), and we 
introduce the corresponding chiral momentum density: 
7r = (f/c 2 . For propagating beam- like fields the integral 
quantities per unit z-length are determined by the 2D 
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spatial integration over the (x, y) plane [9,11]: 
W = J wd 2 r , P = J Pd 2 r , 
X = J X d 2 r , n = J tv d 2 r . 



(19) 



Performing this integration of Eqs. (14)-(17) and using 
/ e lRr d 2 r = (2tt) 2 5 2 (k) along with Eq. (13), we obtain 



W 



r LI J 



<T=±1 



E° 



d 2 k, 



(20) 



meaning of the measure of chirality introduced in [3-5] 
and explains its connection with the polarization state 
and Poynting theorem. Importantly, the chiral momen- 
tum is closely related to the spin angular momentum of 
light, which is represented by the operator S = ak/k 
[9]. Hence, for monochromatic fields the chirality (22) 
and chiral momnetum (23) represent, up to constant 
factors, the averaged helicity of the field and its spin 
angular momentum: 



n = -s . 

c 



(25) 



X 



n 



E 

(7=±1 



E 

a=±l- 



E 

o-=±l" 



E c 



d 2 k 



(21) 





E CT 


d 2 k , 


(22) 
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E CT 


d 2 k . 


(23) 



These equations allow a clear interpretation in terms of 
the quantum-like operators of the corresponding dynam- 
ical quantities [7,9,12]. Introducing the field state vec- 



tor 

( 



= E CT and assuming the convolution defined as 
J2 J d 2 k, we see that Eqs. (20)-(23) can be 

<7=±1 



O 



~E a ) of 



regarded as the expectation values O 

the following operators of energy, momentum, chirality, 
and chiral momentum: 



W = w , P = k, 
X = acj , II = ak . 



(24) 



Thus, using the operator formalism, the 4- 
pscudovcctor of the chirality and chiral momentum, 

^X, 11^ represent just the usual energy-momentum 

(w,k) multiplied by the pseudoscalar of helicity a. This 
remarkably simple result unveils the actual physical 



IV. CONCLUSION 

To summarize, we have examined the measure of the 
optical chirality for a free monochromatic optical field. 
Using the momentum representation (Fourier decompo- 
sition) , we uncovered the close connections of the chiral- 
ity density and chirality flow to the polarization helic- 
ity, energy density, and Poynting energy flow. The ra- 
tio of the chirality density to the energy density reaches 
its maximum absolute value for the fields composed of 
plane waves with the same well-defined helicity. At the 
same time, the local spatial polarization structure result- 
ing from the interference of partial plane waves can be 
rather complicated. Finally, we have determined the in- 
tegral values of the chirality and the corresponding chi- 
ral momentum and have revealed that the usual energy- 
momentum operator multiplied by the helicity underlies 
these quantities. 

Note added.- After submission of this work, a related 
paper [13] came to our attention. Similar quantities and 
conservation laws are considered there, but in the context 
of the Riemann-Silberstein vectors, which enables to deal 
with helicities in the coordinate representation [7b]. 
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